Classical behavior of strongly correlated Fermi systems near a quantum critical point. 

Transport properties 



V. A. Khodel, 1 - 2 J. W. Clark, 2 and M. V. Zverev 1 ' 3 

'Russian Research Centre Kurchatov Institute, Moscow, 123182, Russia 
^McDonnell Center for the Space Sciences & Department of Physics, 

Washington University, St. Louis, MO 63130, USA 
s Moscow Institute of Physics and Technology, Moscow, 123098, Russia 
(Dated: October 22, 2010) 

The low-temperature kinetics of the strongly correlated electron liquid inhabiting a solid is an- 
alyzed. It is demonstrated that a softly damped branch of transverse zero sound emerges when 
several bands cross the Fermi surface simultaneously near a quantum critical point at which the 
density of states diverges. Suppression of the damping of this branch occurs due to a mechanism 
analogous to that affecting the phonon mode in solids at room temperature, giving rise to a classical 
regime of transport at extremely low temperatures in the strongly correlated Fermi system. 
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After a decade of comprehensive studies?^ the preva- 
lence of non- Fermi-liquid (NFL) behavior in strongly cor- 
related Fermi systems is no longer a revelation. How- 
ever, various features of NFL phenomena still await sat- 
isfactory explanation, especially the puzzling observa- 
tions pointing to characteristic classical behavior in a 
quantum-critical regime. For example, at extremely low 
temperatures around 1 mK, experimental measurements 
of the specific heat of two-dimensional (2D) 3 He, as real- 
ized in dense 3 He films are described by the classical 
formula C(T) = j3 + jT, where (3 and 7 are constants. 
Also in contrast to FL theory, the low-temperature resis- 
tivity g(T) of many high-T c compounds in strong mag- 
netic fields, as well as certain heavy- fermion metals, is 
observed^— to vary linearly with T. Hence these systems 
behave as if a major contribution to the collision term 
comes from the electron-phonon interaction, in spite of 
the fact that the phonon Debye temperature To exceeds 
measurement temperatures by a factor 10 2 — 10 3 . 

Previously^ we have attributed the presence of the 
classical term (3 in the specific heat C(T) of 2D liquid 3 He 
to softening of the transverse zero-sound mode (TZSM), 
occurring in the region of the quantum critical point 4 
(QCP) where the density of states N(0), proportional to 
the effective mass M*, diverges. Here we shall address 
the impact of the TZSM on transport properties in the 
QCP regime. The TZSM exists only in those correlated 
Fermi systems where the effective mass M* exceeds the 
bare mass M by a factor more than 3, a requirement 
that is always met on approaching the QCP. In conven- 
tional Fermi liquids, the Fermi surface consists of a single 
sheet, so the TZSM has a single branch with velocity Ct 
exceeding the Fermi velocity vf- Consequently, there is a 
ban against emission and absorption of sound quanta by 
electrons, and the role of the TZSM in kinetics is of little 
interest. However, in heavy-fermion metals, it is usual for 
several bands to cross the Fermi surface simultaneously, 
thereby generating several zero-sound branches. For all 
branches but one the sound velocities are less than the 
largest Fermi velocity. Hence the aforementioned ban is 



lifted, and these branches of the TZSM spectrum expe- 
rience damping, in a situation similar to that for zero- 
spin sound. In the latter instance, Landau damping is 
so strong that the mode cannot propagate through the 
liquid^iiS It will be seen, however, that this is not the 
case for damping of the TZSM, because of the soften- 
ing of this mode close to the QCP. Due to the softening 
effect, the contribution of the damped TZSM to the col- 
lision integral has the same form as the electron-phonon 
interaction at room temperature. On the other hand, we 
will also find that in heavy-fermion metals, analogously 
to the case of liquid 3 He films, softening of the TZSM 
acts to lower a characteristic temperature fit that plays 
the role of a Debye temperature, setting the stage for the 
existence of a classical transport regime at extremely low 
temperatures. 

In the canonical view of quantum phase transitions, 
the QCP has been identified with an end point of the 
line Tjv(p) of a second-order phase transition associated 
with violation of some Pomeranchuk stability condition. 
In turn this violation is associated with divergence of 
the energy derivative <9E(p, e)/de of the self-energy and 
consequent vanishing of the quasiparticle weight z = 
(1 — c?£(p, e)/de)~ 1 in single-particle states at the Fermi 
surface, thus triggerin g 1 ' 2 ' 11 divergence of the effective 
massM* defined by M/M* = z(l + (<9E(p, e)/de°) \ p = pp . 

A number of key experimental studies performed re- 
cently fail to support the canonical view of the QCP. 
In 2D liquid 3 He, experiment^ has not identified any 
phase transition that can be so associated with the point 
of the divergence of N(0). It has been acknowledged^^ 
that a similar situation (cf. Ref. Il3lfl4r ) also prevails for 
the QCPs of heavy-fermion metals. In essence, the point 
where the density of states diverges is separated by an 
intervening NFL phase from points where lines of some 
second-order phase transition terminate. Furthermore, 
these transitions possess unusual properties such as hid- 
den order parameters; therefore within the standard col- 
lective scenario they can hardly qualify as triggers of the 
observed rearrangements. 
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We are therefore compelled to interchange the horse 
and the cart,relative to the canonical scenario. Follow- 
ing Refs. ilalla we attribute the QCP to vanishing of the 
Fermi velocity vp at a critical density p^, which occurs 
if 1 + (<9E(p, e)/9e°) \ P=PF = 0. Accordingly, in this sce- 
nario for the QCP, it is the momentum- dependent part of 
the mass operator that plays the decisive role. 

The authors of most theory articles devoted to the 
physics of the QCP claim that switching on the interac- 
tions between particles never produces a significant mo- 
mentum dependence in the effective interaction function 
/, and hence that the option we propose and develop is ir- 
relevant. This assertion cannot withstand scrutiny. The 
natural measure of the strength of momentum-dependent 
forces in the medium is provided by the dimensionless 
first harmonic F\ = fippM* /tt 2 of the interaction func- 
tion /(pi,P2) of Landau theory. In a system such as 
3D liquid 3 He where the correlations are of moderate 
strength, the result F\ > 6.25 for this measure extracted 
from specific-heat data is already rather large. The data 
on 2D liquid 3 He are yet more damaging to the claim of 
minimal momentum dependence, since the effective mass 
is found to diverge in dense films^^! In the case of QCP 
phenomena occurring in strongly correlated systems of 
ionic crystals, it should be borne in mind that the elec- 
tron effective mass is greatly enhanced due to electron- 
phonon interactions that subserve polaron effects J£r— 

The change in sign of vf at the QCP results not only 
in a divergent density of states, but also in a rearrange- 
ment of the Landau state beyond the QCP. As a rule, 
however, such a rearrangement already occurs before the 
system attains a QCP. This may be understood from 
simple arguments based on the Taylor expansion of the 
group velocity v(p) — de{p)/dp, which has the form 



where vp(pt) = 3v 2 /(8«2), namely 
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Clearly, this transition, identified as a topological phase 
transition 2 ^™ 5 (TPT), takes place already on the disor- 
dered side of the QCP regime where vf(p) is still posi- 
tive. Accordingly, a new hole pocket opens and the Fermi 
surface gains two additional sheets, the new T = quasi- 
particle momentum distribution being given by n(p) = 1 
for p < pi and P2 < P < Pf, and zero elsewhere. The 
emergence of new small pockets of the Fermi surface is an 
integral feature of the QCP phenomenon, irrespective of 
whether the strongly correlated Fermi system is 2D liq- 
uid 3 He, a high-T c superconductor, or a heavy-fermion 
metal. 

At the TPT point pt, the density of states, given by 



N(T) 
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with dv denoting an element of momentum space, is also 
divergent. Inserting the spectrum ([2]), straightforward 
calculation yields N(T — > 0) cx T -1 / 2 , in contrast to the 
behavior N(T — > 0,Poo) °c T~ 2 / 3 obtained^ in the case 
where ui (/>«,) = 0. 

Having tracked the initial evolution of the topology of 
the Fermi surface in the QCP region, our analysis turns 
next to the salient features of the TZSM spectrum in sys- 
tems having a multi-connected (i.e., multi-sheet) Fermi 
surface. We first examine how the TZSM softens in 3D 
systems with a singly-connected Fermi surface, where the 
dispersion relation has the well-known forrn^ 
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in the vicinity of the QCP. We assert that the last coef- 
ficient V2 is positive, to ensure that the spectrum 



e(p) = v f (p)(p~Pf)- 



1 

2p~F~ 



vi(p)(j>-Pf) 



1 

6p F 



2 v 2 (p-Pf) 3 
(2) 

derived from Eq. (flj exhibits proper behavior at large 
separations from the Fermi surface. 

By its definition, the QCP is situated at a density poo 
where vf(p) vanishes. The QCP must in fact correspond 
to an extremum of the function v(p, Poo), which vanishes 
for the first time at p = pp. Thus, the simultaneous van- 
ishing of the coefficient Vx(p 00 ) — (dv(p, Poo)/dp)\ p=PF of 
the second term in the Taylor series is crucial to the oc- 
currence of the QCP. Generally, v\ does not meet this 
additional requirement. However, in relevant cases its fi- 
nite value remains extremely small, making it possible to 
tune the QCP by imposing an external magnetic field. 

When Ui ^ in Eq. ©, this equation unavoidably 
acquires two additional real roots at a critical density pt 



with s = c/v F and F x = fip F M* /tt 2 . The TZSM is seen 
to propagate only if F\ > 6, i.e., M* > 3M. Near the 
QCP where M*(p) — > 00, one has vp/c — > 0, and Eq. flSJ 
simplifies to 
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an analogous formula being obtained for a 2D system. 

To facilitate analysis of damping of the TZSM in sys- 
tems having a multi- connected Fermi surface, we restrict 
consideration to the case of two electron bands. The TPT 
is assumed to occur at one of the bands, so that its Fermi 
velocity, denoted again by vp, tends to zero, while the 
Fermi velocity v a of the other band remains unchanged 
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through the critical density region. The model disper- 
sion relation for the complex sound velocity c = cr + icj 
becomes 
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It can easily be verified that the contribution of the 
second term to the real part of the right-hand side of 
Eq. © is small compared to that of the first term, since 
vf/v toward the QCP. On the other hand, not- 
ing that In [(cr + icj + v )/{cr + icj — v Q )] ~ —in, the 
corresponding contribution ittFiVfCr/ (4uj) to the imag- 
inary part of the right-hand side cannot be ignored, else 
cj = 0. By this reasoning, Eq. ((5J assumes the simplified 
form 
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Importantly, we see then that the ratio ci/cr oc 
y/ M /M*(p) is suppressed in the QCP regime, which al- 
lows us to analyze the contribution of the TZSM to the 
collision term entering the resistivity along the same lines 
as in the familiar case of the electron-phonon interaction. 
By contrast, the group velocities of the damped branches 
of longitudinal zero sound are found to be insensitive to 
variation of the effective mass in the QCP region^iS No 
such quenching by a small parameter M/M* arises, so 
these modes cannot propagate in the Fermi liquid. 

It should now be clear that toward the QCP, the ef- 
fective Debye temperature fl t = w(fc max ) = CRk max goes 
down to zero, independently of the value of the wave num- 
ber fc max characterizing the cutoff of the TZSM spectrum. 
Thus, the necessary condition Q t < T for emergence of 
a regime of classical behavior is always met. However, 
another condition must also be satisfied if there is to ex- 
ist a well-pronounced classical domain at extremely low 
temperature. Consider that the boson contribution 

F B = T fin (l - e- ck/T ) 9(k - k max )dv (11) 

to the free energy is proportional to some power of fc max , 
depending on the dimensionality of the problem. The 
same is true of the corresponding contributions to ki- 
netic phenomena. Therefore the extra condition needed 
is that fc max must not to be too small. We identify fc max 
with the new characteristic momentum arising beyond 
the point of the TPT, namely the distance d = pi — p\ 



between the new sheets of the Fermi surface. Indeed, for 
momenta p situated at distances from the Fermi surface 
significantly in excess of d, the single-particle spectrum 
e(p) is no longer flat. The collective spectrum cj(fc) deter- 
mined from the corresponding Landau kinetic equation is 
no longer soft, and consequently the imaginary part cj of 
c becomes of the same order as cr, barring propagation 
of the TZSM. 

This scenario is illustrated for the key property of re- 
sistivity in electron systems of solids. The kernel of the 
electron-TZSM collision integral underlying the resistiv- 
ity q(T) then contains terms n (p + k)(l — n (p))iV(k) — 
n Q (p)(l - n (p + k))(l + AT(k)) and n (p + k)(l - 
n Q (p))(l + N(k)) - n (p)(l - n D (p + k))JV(k), in which 
iV(k) denotes a nonequilibrium TZSM momentum distri- 
bution and n D (p), a nonequilibrium electron momentum 
distribution of the band that can absorb and emit the 
TZSM. The explicit linearized electron-phonon-like form 
of the corresponding component of the collision integral 
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In this expression, w is the collision probability, No(uj) = 
[exp(o;(k)/T) — 1] _1 is the equilibrium TZSM momentum 
distribution, and Sriij stands for the deviation of the 
real momentum distribution of the electron band labeled 
o from its nonequilibrium counterpart, with = n(jp) 
and rif = n(p + k). In the classical case of the electron- 
phonon interaction, at Td < T one has 8Nq{lo) / du oc 
—T/oj 2 while all the other factors are T- independent, 
resultin g 26 ' 27 in linear variation of the resistivity g(T) 
with T. Based on the analogy we have established be- 
tween the roles of phonons and the TZSM, the resistivity 
of the strongly correlated electron system must obey the 
FL law g(T) oc T 2 only at T < fl t . In the opposite case 
fit < T, the resistivity exhibits a linear dependence on T. 
Imposition of a magnetic field cannot kill the soft mode 
of transverse zero sound as long as the flattening of the 
single-particle spectrum responsible for strong depression 
of the effective Debye temperature fit persists. 

These results and conclusions are in agreement with 
experimental data 6 on the low-temperature resistivity 
of the doped heavy-fermion metal YbRh2(Sio.95Geo.o5)2, 
data which indicates that the linear-in-T dependence of 
the resistivity is robust down to temperatures as low as 
20 mK. Remarkably, the linearity of g(T) continues to 
hold in external magnetic fields up to B ~ 2T, far in 
excess of the critical value B c ~ 0.3T at which this com- 
pound undergoes some phase transition with a hidden or- 
der parameter, 6 A linear T dependence of g(T) is present 
as well in the Hcrtz-Millis spin-density-wave (SDW) sce- 
nario for the QCP in 2D electron systemsj 28 ' 29 However, 
critical spin fluctuations die out at B > B c , since the 
SDW transition is suppressed. Thus the observed behav- 
ior of g(T) contradicts the SDW scenario. Our scenario 
is in fact compatible with this behavior, since the TZSM 
spectrum is less sensitive than the structure of critical 
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spin fluctuations to the magnitude of the magnetic field. 

A linear T dependence of g(T) can also emerge if light 
carriers are scattered by heavy bipolarons . 19 ' 30 However, 
there is no evidence for the presence of these quasiparti- 
cles in heavy-fermion metals. 

Let us now turn to the analysis of the soft TZSM con- 
tribution to the thermopower. One may recall that in the 
classical situation, the phonon-drag thermopower Sd(T) 
associated with nonequilibrium phonons is know n 26 i 27 i 31 
to account for a substantial part of the Seebeck coefficient 
S(T). The same is true for the class of quantum-critical 
systems considered here, except that the domain in which 
the drag term Sd(T) contributes appreciably extends 
down to extremely low temperatures. Significantly, at 
T — > the drag contribution increases as T 3 , whereas at 
T > fi t it falls off as T _1 , producing a bell-like shap o 26 i 27 
of Sd(T) with a sharp maximum at T ~ fl t . Since the re- 
maining contributions to S arc rather smooth, this salient 
feature of Sd appears to be responsible for the change 
of sign of the full Seebeck coefficient S(T) at extremely 
low T observed experimentally^ 2 - in the heavy-fermion 
metal YbRli2Si2, as well as the nonregular behavior of 
S(T 0) found in several heavy-fermion compounds. 33 



The TZSM scenario proposed here predicts that the See- 
beck coefficient in YbRh2(Sio. 95660.05)2 will exhibit the 
same anomalous behavior at magnetic fields substantially 
exceeding a corresponding critical value B c . 

In summary, we have investigated the conditions pro- 
moting the formation of soft damped collective modes 
that play the same role in kinetic phenomena as phonons. 
We have shown that such a damped soft branch belonging 
to the transverse zero-sound mode emerges when several 
bands cross the Fermi surface simultaneously, with one of 
the bands subject to a divergence of the effective mass of 
carriers. We have discussed prerequisites for lowering the 
corresponding Debye temperature fit and demonstrated 
that the inequality tt t < T is met near the quantum crit- 
ical point, where a classical regime sets in at extremely 
low temperatures. 
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